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While standard scaling arguments show that a system of non-interacting electrons in two dimen- 
sions and in the presence of uncorrelated disorder is insulating, in this work we discuss the case 
where inter-impurity correlations are included. We find that for point-like impurities and an infi- 
nite inter-impurity correlation length a mobility edge exists in 2D even if the individual impurity 
potentials are random. In the uncorrelated system we recover the scaling results, while in the in- 
termediate regime for length scales comparable to the correlation length, the system behaves like 
a metal but with increasing fluctuations, before strong localization eventually takes over for length 
scales much larger than the correlation length. In the intermediate regime, the relevant length 
scale is not given by the elastic scattering length but by the inter-impurity correlation length, with 
important consequences for high mobility systems. 



It is generally believed that a non-interacting two di- 
mensional (2D) system in the presence of disorder is al- 
ways insulating. This result is based on extensive work 
on scaling theory pioneered by the 'gang of four' It 
has also received considerable amount of numerical sup- 
port |2|. The general statement can be summarized as 
follows. In scaling theory the localization length in 2D is 
given by L c ~ \ e Kk F \/2 jgj ^ w j iere \ j s \\^ e elastic scatter- 
ing length. Hence, as soon as A is finite, the localization 
length is finite. However, because A appears in the ex- 
ponent, when fcfA >> 1 the localization length can be 
extremely large and difficult to probe experimentally. In 
the standard Born approximation A -1 ~ nj, where ni 
is the two-dimensional impurity concentration, hence, a 
non-zero ni would lead to a finite L c . At zero tempera- 
ture a finite localization length implies that the resistance 
diverges exponentially for a system size exceeding the lo- 
calization length. Equivalently, for an infinite system, 
the resistance diverges when the temperature, T tends 
to zero because the phase coherence length 1$ is infinite 
at zero T. 

When considering discrete models based on Anderson's 
disordered tight binding model Q very similar results 
are obtained. All states are localized for any strength 
of disorder. These results apply only if all sites are un- 
correlated. Indeed, even in ID there exist special long 
and short range correlations in the disorder, which can 
lead to the existence of extended states in these systems 
3- Similarly special systems can also be found in 2D 
6j. The exact conditions under which these localization 
conditions apply has recently gained considerable inter- 
est because 2D electronic systems confined in a variety 
of semiconducting structures, such as Si-MOSFETs and 
p-type GaAs/AlGaAs have shown a strong metallic-like 
T-dependence of the resistance 0, @ down to the low- 
est experimental T. In metallic, we understand a positive 
or vanishing derivative of the resistivity as a function of 
the temperature, i.e., dp/dT > 0. 

Because in these systems disorder is always present, 



which is expected to lead to localization, most of the re- 
cent theories have considered interactions between elec- 
trons 10;] as a possible mechanism for the metallic be- 
havior. In this work, we consider instead the case of cor- 
relations between impurities, which leads to a genuine 
metal-insulator transition for an infinite inter-impurity 
correlation length I — * oo for the non-interacting 2D dis- 
ordered system (with random potentials). Moreover, for 
a finite I, the crossover to localization is determined by I 
and not A. 

We start by considering the standard disordered tight 
binding Anderson model in 2D: 
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When V n . m = we recover the free particle case with 
solutions of the type ip n , m = e ikn e tpm and eigenvalues 
E = 2cos(fc) + 2cos(p). When V n>m is random and un- 
correlated, scaling theory predicts that all states are lo- 
calized for any disorder strength. The same result holds 
in ID. However, when the potential is correlated the sit- 
uation becomes very different. In ID, when only every 
2nd site is random, i.e., V% n — and Vm+\ is random 
ipn = cos(n7r/2) is an extended solution of the disordered 
system, with energy E = 0. In this model only one en- 
ergy has an infinite localization length and since at all 
other energies all states are localized there is no band of 
extended states jjjj- However, when we consider a simi- 
lar potential in 2D something very interesting happens. 

When defining Vin,m — and V2 n +i, m as random, we 
can write a solution to equ. 1 as ip n ,m = cos(n7r/2)e lfcm , 
corresponding to E — 2cos(£;). Clearly, this is a fully 
extended state. The difference in the 2D case compared 
to ID is that we now have a band of extended states 
for — 2 < E < 2. For every energy in this band there 
is exactly one wave- function, which can be written in 
this form, hence there is no degeneracy. Therefore, when 
evaluating the two-terminal conductance of the system 



2 




one obtains G = 2e 2 /h for -2 < E F < 2 and G = 
otherwise, since all other states are localized. This gives 
rise to a band of extended states around the band center 
very similar to the 3D Anderson model, which leads to a 
metal-insulator transition when the Fermi energy crosses 
E = —2 and E = 2. Similarly, when only the V an ^mS 
are random, but all other potentials are zero, the system 
has conductance steps at E — ±2 + 2cos(?7r/a), with 

1 = 1, 2, 3, . . . , a — 1 and a an integer a > 2. We have 
evaluated the conductance numerically as a function of 
the Fermi energy and the results are presented in fig. 1. 
In order to calculate the temperature dependence of the 
resistance we simply used G(T) — J dEf T (E — Ep)G(E), 
where fx is the Fermi-Dirac distribution function. For 
non-interacting electrons this is the dominant tempera- 
ture dependence since inelastic scattering is strongly sup- 
pressed at low temperatures. Moreover, because Ep is 
equivalent to the density, this system exhibits a genuine 
metal-insulator transition as a function of density. 

The continuous case can be treated in an analogous 
way by considering <5-impurities with random amplitudes 
on a lattice, which is similar to the ID case discussed in 
ref. 0. In this case the random potential is V(x,y) = 
Y^ n , m Vn,m(y)S(na-x)6(ma-y), where V n , m are random. 
The special energies for the conductance steps are now 
given by E = n 2 in units of 7i 2 7r 2 /2ma 2 , where a is the 
corresponding lattice constant of the random impurities 
and n — 1,2, • • •. The corresponding conducting wave- 
functions (extended states) are simply given by ip(x, y) — 
sin(rnrx/a)e lky for E = n 2 + (ka/n) 2 . In the inset of fig. 

2 the zero temperature conductance, G, is shown as a 
function of the Fermi energy. 

The temperature dependence is shown for different val- 
ues of Ep in fig. 2. Hence, in this simple model in which 
6 impurities are on a lattice with random amplitudes, 
we have a Ep induced metal-insulator transition at Eq, 
independent of the disorder strength. The temperature 
dependence was obtained by using the Fermi-Dirac distri- 
bution R- 1 = J dEf T (E- Ep)G(E). It is quite striking 
to note that the overall shapes of the curves are very sim- 



FIG. 2: Temperature dependence of the 2-terminal resistance 
for various values of Ef- The inset shows the zero tempera- 
ture dependence of the conductance as a function of Ef and 
the dotted line illustrates the overall square root dependence 
of the conductance steps 

ilar to the experimental ones in high mobility 2D systems 
|8|. The temperature scale in these experiments is typi- 
cally close to the Fermi temperature Tp, which is of the 
order of IK and very similar to the relevant scale in our 
simple model. 

In the above discussion we considered the two-terminal 
conductance of the system. Experimentally, this is the 
quantity which is normally measured. However, in order 
to neglect the contact resistance, four terminal measure- 
ments are typically performed, but at zero magnetic field 
this is equivalent to a two terminal measurement assum- 
ing zero contact resistance. In addition, in 2D the value 
of the conductance is converted to conductivity by tak- 
ing into account the geometry of the sample. Hence the 
conductivity in our case would be a = G/n, where □ is 
the number of squares or geometrical factor. Since the 
conductance of our model does not depend on the geom- 
etry or length of the sample, our system would have an 
infinite conductivity if the sample was infinitely long but 
finite in width and of zero conductivity if infinitely wide 
and of finite length. At first this does not seem physical 
but if one considers that each experimental system has a 
finite I,), at non-zero T, which is equal in every direction 
for an isotropic system, the relevant geometry is close to a 
square. This is very reasonable because for any practical 
comparison with experiments, quantum effects such as 
localization can only occur within 1$. Indeed, we could 
model the experimental system as a network of quan- 
tum coherent squares coupled classically and we would 
recover a — G, which is the expected result in which the 
value of the critical conductivity does not depend on the 
geometry. This normalization was used for all figures. 

Experimentally, a large applied parallel magnetic field 
tends to exponentially increase the resistivity |l3| . In a 
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non-interacting picture a parallel field would simply shift 
the two spin subbands and therefore not significantly al- 
ter the overall behavior, although it could decrease the 
critical transition density. However, when comparing to 
experimental systems in a parallel field other effects be- 
come very important too. For instance, the 2D plane is 
not entirely flat because there are micron sized in-plane 
deviations, which are of the order of 4 to 30nm high ver- 
tically ^3 > which can allow random fluxes to penetrate 
the 2D system and induces additional localization. These 
effects are beyond the scope of this work and will not be 
considered further. 

Coming back to our model, where the disorder is con- 
fined to a regular lattice, we could show that there is 
a well defined metal-insulator transition, very similar to 
the experimental situation. However, impurities do not 
usually form a perfect crystal during the growth process. 
They are more likely to be somewhat correlated but not 
to the extent to form a lattice. In order to generalize 
our model, we start by considering the case, where the 
impurities are uncorrelated I = and where the density 
of impurities is low before we turn to the more general 
case of arbitrary I. 

For the numerical calculation we considered the tight- 
binding equation (1) with finite extent L in the n direc- 
tion and infinite in m. Without disorder we have plane 
wave solutions of the type ipn,m = sm(npir / L)e lkm , with 
E = 2 cos(p7r/L)+2 cos(fc). This system can be viewed as 
a quasi-lD quantum wire, where the degeneracy of states 
depends on L and E and the conductance is simply 2e 2 /h 
multiplied by the degeneracy. In the middle of this sys- 
tem we now add Nj punctual impurities in a rectangle of 
size (width=L)x (length=L). Since contacts are made by 
diffusing a metal into the 2D layer, the metal side of the 
contact has a much larger electron density than the 2D. 
Therefore, to model the effect of contacts we have consid- 
ered the following situation. Instead of using a source and 
drain region of the same size as the disordered region, we 
consider a much larger (essentially infinite) contact point 
to the disordered region. The model is illustrated in fig. 
3, where the source and drain region has no disorder and 
is much wider than the disordered region. Hence, the 
wide regions correspond to the metal contacts, whereas 
the disordered narrow region corresponds to the 2D sys- 
tem under study. We evaluated numerically the conduc- 
tance of the system using standard transfer matrix tech- 
niques and show the conductance as a function of the 
number of impurities for different sizes of the disordered 
region in fig. 3. In our model we used a square geome- 
try L x L, which has the particularity that at the band 
center the conductance G is independent of L for V = 
and a = G. The curves were obtained by averaging over 
ten different disorder configurations. 

Fig. 3, clearly shows that a only depends on the num- 
ber of impurities and not on the size of the system. Hence 
the dependence is simply <r(iVj) = a(L nj). This depen- 
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FIG. 3: The conductivity as a function of the number of im- 
purities for different sizes and values of the impurity potential, 
with -0.5 < V x < 0.5, -1< V 2 < 1 and V 3 = 1. 




FIG. 4: The conductivity as a function of the number of im- 
purities for different sizes and for — 1 < Vi < 1. Here only 
every second site is disordered. 

dence holds for the entire range until a saturates to the 
value where ni ~ 1. At large enough number of impuri- 
ties and in the strong localization regime we have an ex- 
cellent fit to the expression a(Ni) ~ &- a ^ ^ e -a^njL ^ 
where a depends only on the impurity strength and 
L^ 1 = a^/nj. The quality of the fit and the simplicity 
of the fitting expression is quite remarkable. This result 
also implies that for any non-zero uncorrelated impurity 
density the system is always insulating, which is consis- 
tent with scaling analysis. 

However, when the impurities are distributed on a lat- 
tice, i.e., I = oo, like in fig. 1, we also obtain a universal 
dependence only on Nj independent of L, i.e., a(Ni). 
The result is shown in fig. 4. The only difference here 
is that a(Nj — ► oo) = 2e 2 /h instead of zero. For this 
data set we considered that only every second site is ran- 
dom, which in this case leads to a metallic conductance 
of 2e 2 /h at the band center. This data set was obtained 
without averaging over different disorder configurations 
and each data point represents an independent impurity 
configuration. 
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FIG. 5: The conductivity as a function of the number of im- 
purities for different ratios of correlated to uncorrelated im- 
purities. The straight line is the fit a ~ exp{—aN^ 2 ). 

We now consider the most interesting case, where we 
change the inter- impurity correlation length I. The re- 
sults in fig. 5 illustrate three cases (I = Q, I ^ 0, 1 = oo) 
for a fixed system size L since a does not depend explic- 
itly on L but only on Nj. In order to vary I we chose to 
introduce two types of impurities Nj orr and N™ ncorr and 
Nj = N c j orr + Nf ncorr . The Nf orr impurities will fall 
randomly on the lattice (of spacing 4 for fig. 5), whereas 
the N™ ncorr impurities are not correlated to any site. 
The correlation length is therefore given by the density 
of uncorrelated impurities I = £ 2 /iV" ncor ''. 

Overall, the conductivity remains metallic for length 
scales comparable to the correlation length and the con- 
ductance fluctuations are close to zero in this metallic 
regime. For length scales exceeding the correlation length 
we first have increasing conductance fluctuations before 
the system eventually localizes at large enough length 
scales L ^S> I with a localization length equal to the uncor- 
related case. The fluctuations in the intermediate regime 
are much larger than for the uncorrelated system and 
constitute an important signature of the correlation in 
the "metallic" -like state. Interestingly, a very similar en- 
hancement of the conductance fluctuations is seen close 
to the transition from the metallic to insulating behav- 
ior in Si-MOSFET systems and GaAs [II]. In our model 
of correlated impurities it is therefore the inter-impurity 
correlation length I, which controls the range of metal- 
lic behavior and not the semiclassical elastic scattering 
length A, which can be much smaller than I. Hence, a 
large I can significantly enhance the range of metallic 
behavior in 2D, which constitutes one of the the main re- 
sults of this article. This is in stark contrast to the case 
where the impurity range is increased. Indeed, in this 
case it was shown that a larger impurity range enhances 
localization if impurities remain uncorrelated [l6|. The 
situation with correlated impurities was not considered. 



Hence, this shows that inter-impurity correlations sup- 
press localization whereas large impurity ranges enhance 
localization. In an applied perpendicular magnetic field, 
it was recently shown that a long-range disorder poten- 
tial might induce a metal-insulator transition |17| too. 
Clearly, correlations between impurities have a signifi- 
cant impact on the localization properties and have to 
be considered when comparing to experimental systems. 

In conclusion, we have analyzed the impact of inter- 
impurity correlations on the localization properties of a 
non-interacting 2D electronic system. For impurities on 
a lattice and with random potentials, corresponding to 
an infinite inter-impurity correlation length, the system 
exhibits a true metal-insulator transition. A finite cor- 
relation length enhances the scale over which the system 
is conducting before the system is eventually localized. 
A large correlation length can explain the large conduc- 
tance fluctuations and metallic behavior of experimental 
systems. We also showed that the conductivity is only 
a function of the number of impurities, Ni and that the 

— 1/2 

localization length is given by L c ~ Nj for L^> I. 
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